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Abstract 

In  addition  to  the  efforts  to  reduce  the  costs  of  renewable  energy  technologies  and  electrolyzers,  the  development  of  suitable  controllers 
are  needed  for  cost-competitive  electricity  production  by  renewable-hydrogen  power  plants.  In  this  paper,  a  novel  control  is  proposed  for  a 
wind-electrolysis  system,  which  match  the  wind  power  output  to  the  electrolyzer  power  requirements,  thus  gaining  in  system  performance. 
It  basically  consists  in  continuously  shaping  the  power  reference  of  a  conventional  maximum  power  point  tracking  algorithm.  Thus,  high 
aerodynamic  power  conversion  efficiency  is  achieved  fulfilling  at  the  same  time  the  electrolyzer  specifications.  This  control  strategy  is 
developed  using  concepts  of  the  reference  conditioning  technique  and  of  the  sliding  mode  control  theory.  The  proposed  control  algorithm  is 
extremely  simple,  easy  to  implement  and  robust  to  parameter  uncertainties. 
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1.  Introduction 

Clean  and  sustainable  electricity  generation  from  renew¬ 
able  resources  have  experienced  a  significant  development 
during  the  last  decades.  Nowadays,  alternative  energy  conver¬ 
sion  systems  are  commonly  found  in  different  autonomous 
and  grid-connected  applications.  However,  the  intermittent 
and  seasonal  variability  of  the  primary  resources  hampers  a 
greater  penetration  of  these  clean  energies  into  the  power  gen¬ 
eration  market.  The  aggregation  of  different  types  of  energy 
resources,  for  instance  wind  and  solar,  helps  to  attenuate  this 
limitation.  Nevertheless,  the  best  solution  consists  in  storing 
the  energy  captured  in  excess  to  be  used  during  short-supply 
periods.  Although  batteries  are  extensively  used  in  low  power 
systems,  their  use  in  medium  and  large  power  applications  is 
inviable  because  of  their  relatively  low  storage  density  (en¬ 
ergy/volume).  A  feasible  option  is  the  production  of  hydro¬ 
gen  through  water  electrolysis  for  its  subsequent  use  either  to 
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produce  electricity  or  to  supply  fuel  engines.  In  fact,  hydro¬ 
gen  allows  both  storage  and  transportation  of  large  amounts 
of  energy  at  much  higher  energy  densities  [1-3].  Traditional 
storage  options  are  actually  viable  for  hydrogen  consump¬ 
tion  in  large  plants.  Yet  the  on-board  storage  for  use  on  light 
vehicles  is  an  open  problem  being  actively  investigated  [4]. 

Due  to  the  comparatively  low  cost  of  wind  technology, 
along  with  the  persistent  growth  in  installed  wind  power 
capacity  over  the  world,  wind-electrolysis  is  the  favorite 
candidate  to  become  the  first  economically  viable  renewable- 
hydrogen  production  system.  Before  this  can  be  realized, 
however,  optimization  of  the  whole  system  is  needed 
in  addition  to  the  ongoing  efforts  to  reduce  the  costs  of 
electrolysis  technology.  Many  research  projects  and  articles 
reported  in  the  literature  about  hydrogen  production  plants 
are  concerned  with  the  sizing,  economics  and  power  flow 
management  of  the  devices  comprising  the  system  (wind 
turbine,  electrolyzer,  batteries,  etc.),  in  order  to  guarantee 
service  quality  at  minimum  cost.  These  problems  are  closely 
associated  to  the  local  resource  characteristics  [5-8]. 

Two  different  approaches  are  found  with  regards  to  the 
configuration  of  stand-alone  wind-electrolysis  system.  In 
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most  cases,  the  turbine  and  electrolyzer,  with  their  own  ded¬ 
icated  power  electronics  and  controllers,  are  connected  to 
a  constant-voltage  dc-bus.  Such  a  decoupled  configuration 
presents  the  advantage  that  the  power  converters  of  the  tur¬ 
bine  and  the  electrolyzer  can  be  controlled  separately.  In 
this  scheme,  the  variable-frequency  variable-voltage  ac  at  the 
wind-driven  generator  terminals  is  converted  to  the  constant 
dc-bus  voltage  through  an  ac-dc  converter,  whereas  a  dc-dc 
converter  takes  the  dc-bus  voltage  and  provides  a  dc  voltage 
suitable  for  electrolyzer  operation.  The  voltage  of  this  de¬ 
bus  is  fixed  by  a  battery  pack,  which  also  has  to  manage  the 
power  mismatches  between  the  wind  turbine  and  the  elec¬ 
trolyzer  and  serves  as  a  filter  of  the  input  power  fluctuations. 

A  more  recent  approach  proposes  eliminating  the  compo¬ 
nent  duplication  in  the  interface  between  the  turbine  and  the 
electrolyzer.  In  fact,  the  pair  of  power  converters  and  the  de¬ 
bus  can  be  replaced  with  a  single  ac-dc  converter  taking  the 
ac  generator  voltage  and  providing  an  appropriate  voltage 
to  the  electrolyzer.  Thus,  the  efficiency  of  the  energy  con¬ 
version  can  be  increased  and  the  overall  cost  of  the  system 
can  be  appreciably  reduced  [9,10].  Obviously,  this  coupled 
configuration  requires  the  development  of  controllers  specif¬ 
ically  designed  for  this  application.  In  this  context,  the  cur¬ 
rent  paper  addresses  the  control  of  a  wind-powered  hydrogen 
production  system  to  match  the  input  wind  power  to  the  elec¬ 
trolyzer  requirements,  thereby  increasing  the  captured  energy 
and  improving  the  dynamic  performance.  Thus,  the  sizing  of 
the  system  required  to  accomplish  the  specifications  can  be 
substantially  reduced. 


2.  System  description 

2.1.  The  plant 


The  system  under  consideration  is  sketched  in  Fig.  1.  Its 
main  components  are: 


•  Wind  turbine.  The  dynamic  behavior  of  the  wind  turbine 
can  be  described  by  the  first-order  differential  equation: 


d  T2 
dr 


'j(Tj  -  7b) 


(1) 


where  Q  is  the  shaft  speed,  J  the  inertia  of  the  rotating 
parts,  Tq  the  electric  torque  developed  by  the  generator 
and  Tj  is  the  aerodynamic  torque: 


product  of  this  latter  coefficient  and  the  wind  power  [11]: 
/MA,  w )  =  -pjrr2Cp(X)w2.  (3) 

As  it  can  be  observed  in  Fig.  2,  the  power  conver¬ 
sion  efficiency  takes  its  maximum  Cpopt  at  the  optimum 
tip-speed-ratio  kopt.  So,  fixed-speed  wind  turbines  oper¬ 
ate  with  maximum  efficiency  at  just  one  wind  velocity. 
Therefore,  variable-speed  controllers  are  commonly  im¬ 
plemented  in  order  to  track  the  optimum  power  Popt(w)  = 
TM^ojib  w).  This  sort  of  control  strategy  is  called  maxi¬ 
mum  power  point  tracking  (MPPT)  controller. 

•  Permanent  magnet  synchronous  generator.  For  the  pur¬ 
poses  of  this  work,  the  generator  can  be  modeled  as  a 
three-phase  sinusoidal  voltage  source  EF  in  series  with  a 
synchronous  reactance  As  as  shows  the  per-phase  circuit 
sketched  in  Fig.  3.  The  phase  current  and  voltage  at  the 
generator  terminals  are  7s  and  Vs,  respectively.  Note  that 
the  terminal  voltage  Vs  slightly  differs  in  magnitude  and 
phase  from  the  emf  EF  because  of  the  voltage  drop  across 
As  when  the  current  7s  passes  through  it.  The  mechanical 
variables  of  the  generator,  i.e.  speed  and  torque,  are  related 
to  the  electrical  ones  as  follows: 

Ef  =  -0T2 

2 

3  Is  Vs  cos  tp  1-1  ’ 

T°  = - a — 

where  p  is  the  number  of  poles  of  the  machine,  0  the 
stator  flux  linkage  and  <p  is  the  phase  angle  between  Vs  and 
Is-  The  former  equation  is  derived  from  electromagnetic 
considerations,  whereas  the  latter  is  obtained  from  power 
balance  between  mechanical  ( Pq  =  T(,Q)  and  electrical 
(Ps  —  37s  Vs  cos  0)  powers  [12]. 

•  Electronic  converter.  The  electronic  devices  are  aimed  at 
converting  the  ac  three-phase  speed-dependent  generator 
voltage  Vs  into  a  dc  voltage  Vh  suitable  for  hydrogen  pro¬ 
duction.  As  it  was  previously  mentioned,  a  single  ac-dc 
converter  is  considered  in  this  paper  as  interface  between 
the  turbine  and  the  electrolyzer.  That  is,  the  voltage  at 
the  generator  terminals  is  directly  converted  to  the  voltage 
at  the  electrolyzer  terminals.  The  output  and  input  volt¬ 
ages  of  the  converter  are  related  by  the  controlled  factor 
<5,  which  is  associated  to  the  duty  cycle  of  the  electronic 
switching  devices.  Since  incoming  and  outgoing  powers 
are  matched,  i.e.  P$  =  P]\,  input  and  output  currents  have 
a  similar  ratio.  That  is, 


Tt(X,  w)  = 


1 

2 


pitr2  Ct(X)w2  . 


(2) 


5VS  =  VH 
<57h  =  37s  cos  <p. 


(5) 


This  aerodynamic  torque  is  a  nonlinear  function  of  the 
wind  velocity  w  and  the  tip-speed-ratio  X  =  rT2/w,  with 
r  being  the  blades  length  and  p  being  the  air  density.  Fig. 
2  displays  the  torque  coefficient  Ct  (dashed  line)  as  well 
as  the  power  efficiency  Cp  =  /.C'x  (solid  line).  Hence,  the 
power  captured  by  the  turbine  can  be  expressed  as  the 


It  is  assumed  in  this  paper  that  cos  (p  =  1 ,  i.e.  that  the  input 
voltage  and  current  are  in  phase.  This  is  the  case  in  many 
ac-dc  converters,  such  as  those  having  an  uncontrolled 
rectifier  at  the  input  [12]. 

•  Electrolyzer.  The  electrolyzer  is  modeled  by  an  empir¬ 
ically  obtained  current-voltage  (7h  —  Vh)  curve.  This 
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Fig.  1 .  Wind-powered  hydrogen  plant  scheme. 


Fig.  2.  Torque  (dashed)  and  power  (solid)  coefficients  of  the  turbine. 

In  —  Vn  characteristic  is  influenced  by  the  electrolyte  tem¬ 
perature  T  and  can  be  approximated  by  a  logarithmic  law 

[13]: 

Vn  =  Vi  +  n(T)hi  +  m(T ) lg  ^1  +  /r)  (6) 

where  Vr  is  the  reversible  cell  voltage  and  n(T),  m{T)  and 
q(T )  are  quadratic  functions  of  T.  Fig.  4  plots  the  /h  —  Vh 
curves  of  a  2.25  kW  von  Hoerner  electrolyzer  for  two  limit 
temperatures  [14]. 


Fig.  4.  Current-voltage  characteristic  of  the  electrolyzer. 


trolyzer.  In  fact,  this  subsystem  requires  a  special  treatment 
because  it  is  the  most  affected  by  wind  resource  variability. 
The  main  control  objective  is  the  maximization  of  the  energy 
captured  by  the  turbine  satisfying  at  the  same  time  the  fol¬ 
lowing  requirements  for  proper  operation  of  the  electrolyzer 
[15]: 

( 1 )  Current  and  voltage  must  be  upper-bounded  at  their  rated 
values,  /N  and  Vr,  respectively.  In  fact,  though  operation 
at  high  current  levels  is  desired  to  improve  the  efficiency 
of  the  electrolyzer,  the  system  must  be  maintained  within 
its  safe  operation  area. 

(2)  A  minimum  current  Im  must  be  guaranteed  with  minimal 
connection  to  the  auxiliary  power  supply  because  oper¬ 
ation  at  low  currents  affects  the  quality  of  the  product. 

(3)  The  wind  power  must  be  smoothed  before  being  supplied 
to  the  electrolyzer.  As  reported  in  some  recent  studies, 
variability  of  the  electrolyzer  power  may  increase  the  in¬ 
ternal  wear  as  well  as  the  impurities  and  energy  losses. 
So,  limiting  the  rate  of  change  of  the  input  power  is  de¬ 
sired. 


2.2.  The  controller 

Stand-alone  hydrogen  production  systems  need  several 
control  loops  as  well  as  a  supervisory  control.  In  this  pa¬ 
per,  we  focus  on  the  control  of  the  subsystem  comprising  the 
wind  turbine,  the  generator,  the  ac-dc  converter  and  the  elec- 


Fig.  3.  Per-phase  model  of  the  permanent  magnet  synchronous  generator. 


Note:  Although  these  items  are  established  to  suitably  sup¬ 
ply  an  electrolyzer  with  wind  power,  similar  specifications 
could  be  additionally  stated  for  the  battery  pack.  Moreover, 
though  the  proposed  control  strategy  is  particulary  focused 
on  wind-powered  electrolyzers,  it  is  also  applicable  to  other 
renewable-electrolysis  systems. 

Fig.  5  shows  a  schematic  diagram  of  the  control  system 
proposed  in  this  paper.  It  basically  consists  of  an  MPPT  con¬ 
troller  conventionally  used  in  wind  turbine  control  and  a  novel 
power  reference  conditioning  loop  proposed  here  to  accom¬ 
plish  the  specifications  of  the  electrolyzer.  The  proposed  con¬ 
ditioning  circuit  is  enclosed  with  dashed  line  in  Fig.  5. 

Let  ignore  for  the  moment  the  conditioning  circuit  and 
focus  on  the  MPPT.  That  is,  let  suppose  that  Pr  =  P*.  The 
MPPT  comprises  an  outer  control  loop  and  an  inner  control 
loop.  The  former  computes  the  optimum  power  reference 
as  function  of  the  shaft  speed,  assuming  that  the  turbine  is 
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Fig.  5.  Scheme  of  the  control  system  with  the  proposed  reference  condition¬ 
ing. 

operating  at  the  optimum  tip-speed-ratio: 

P\d)  =  P0pt(/-.aA0pt)  =  \pnr5C^£2\  (7) 

1  ^opt 


age  and  current  magnitudes.  Hence,  the  captured  energy  is 
increased  while  the  requirements  of  the  electrolyzer  are  sat¬ 
isfied.  Reference  conditioning  is  carried  out  in  a  controlled 
fashion,  thus  preserving  stability.  This  proposal  for  the  con¬ 
trol  of  the  wind-powered  electrolyzer  is  the  main  contribution 
of  the  paper  and  is  described  in  more  detail  in  the  following 
section. 


3.  Electrolyzer  power  conditioning  control  strategy 

In  this  section,  a  control  algorithm  to  satisfy  the  additional 
specifications  enumerated  above  for  the  electrolyzer  is  pre¬ 
sented.  It  consists  in  shaping  the  power  reference  via  a  dis¬ 
continuous  control  action  using  tools  of  sliding  mode  control 
(SMC)  theory.  Reference  conditioning  using  SMC  concepts 
was  originally  proposed  to  design  anti-windup  techniques  for 
constrained  process  control  [16].  In  this  case,  similar  ideas 
are  used  to  limit  the  operating  ranges  of  the  electrolyzer. 


The  latter,  i.e.  the  inner  loop,  adjusts  the  duty  cycle  of  the 
electronic  converter  in  order  to  track  the  reference  power  Pr. 
In  this  loop,  a  conventional  PI  controller  is  usually  employed 
to  cancel  steady  state  errors.  As  a  result,  the  turbine  speed  £ 2 
converges  to  the  optimum  value  £2opt(u;)  =  A.optw/r.  Hence, 
the  maximum  power  point  (^opt(w),  P0pt(w))  is  successfully 
tracked  with  the  closed-loop  dynamic  response  dominated  by 
the  drive-train. 

MPPT  control,  with  its  inherent  variable-speed  operation, 
has  an  additional  advantage  with  respect  to  fixed-speed  op¬ 
eration.  In  fact,  the  more  flexible  connection  to  the  electric 
grid  provides  damping  at  the  resonance  modes  of  the  drive- 
train,  thus  reducing  high  frequency  oscillations  of  torque  and 
power.  Despite  this  power  quality  improvement,  MPPT  con¬ 
trol  strategies  usually  violate  the  particular  specifications  of 
the  electrolyzer.  In  fact,  the  captured  power  may  exceed  the 
rated  power  of  the  electrolyzer  and  must  be  limited.  Addi¬ 
tionally,  the  torque  efforts  developed  to  track  the  optimum 
power  point  translate  into  fast  changes  of  the  electric  power 
that  are  transmitted  to  the  electrolyzer  through  the  fast  elec¬ 
tronic  devices.  This  is  not  an  easy  problem  to  solve.  On  one 
hand,  smoothing  the  power  reference  P*(£2)  via  a  low-pass 
filter  may  lead  to  undesirable  transient  responses.  Moreover, 
the  turbine  may  slow  to  a  stop  when  the  wind  velocity  drops. 
On  the  other  hand,  slowing  down  the  response  of  the  inner 
loop  may  make  the  overall  system  unstable  too.  Then,  a  fea¬ 
sible  option  to  overcome  this  problem  would  be  the  design 
of  an  MPPT  controller  taking  the  electrolyzer  specifications 
into  account  from  the  beginning.  This  is  a  difficult  task  be¬ 
cause  of  the  nonlinear  behavior  of  all  the  devices  involved  in 
the  system. 

In  this  paper,  a  completely  different  approach  is  explored. 
In  fact,  ideas  of  power  tracking  control  are  still  used,  but  now, 
an  auxiliary  reference  conditioning  loop  is  incorporated.  The 
aim  of  this  loop  is  shaping  the  power  reference  in  order  to 
limit  the  rate  of  change  of  the  input  power  as  well  as  the  volt¬ 


3.1.  Background  on  sliding  mode  control 

A  variable  structure  system  comprises  a  set  of  continuous 
subsystems  with  an  associated  switching  function  that  de¬ 
termines  a  manifold  on  the  state  space,  the  so-called  sliding 
surface.  According  to  the  sign  of  the  switching  function,  the 
control  signal  takes  one  of  different  possible  values,  leading 
to  a  discontinuous  control  law.  The  basic  idea  is  to  enforce  the 
state  to  reach  the  prescribed  sliding  surface  and,  henceforth, 
to  slide  on  it  through  a  very  fast  switching  action.  Once  this 
particular  mode  of  operation  is  established,  known  as  sliding 
mode  (SM)  or  sliding  regime,  the  prescribed  manifold  im¬ 
poses  the  new  system  dynamics.  Among  other  attractive  fea¬ 
tures,  sliding  regimes  are  easy  to  implement,  reduce  the  order 
of  the  system  dynamics,  and  provide  robustness  to  matched 
uncertainties  and  external  disturbances  [17-19]. 

Consider  the  following  dynamical  system: 

77  =  fix)  +  g(x)u  (8) 

at 

where  x  e  Sli"  is  the  system  state,  u  the  discontinuous  control 
signal,  and  f(x )  and  g(x)  are  vector  fields  in  ilf" .  The  variable 
structure  control  law  is  defined  as: 

u~.  if  a(x)  <  0 

’  (9) 

m+,  if  cr(.r)  >  0 

according  to  the  sign  of  the  auxiliary  output  er(.r).  The  slid¬ 
ing  surface  S  is  defined  as  the  manifold,  where  the  auxiliary 
output,  also  called  switching  function,  vanishes.  That  is, 

5={rer|ffW  =  0}.  (10) 

When,  as  a  result  of  the  switching  policy  (9),  the  reaching 
condition: 

lim  aa  <  0 

cr— >0 


(ii) 
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locally  holds  at  both  sides  of  the  surface,  a  switching  sequence 
at  very  high  frequency  (ideally  infinite)  occurs,  constraining 
the  system  state  trajectory  to  slide  on  S. 

For  a  sliding  motion  to  exist  on  S  (in  other  words,  to  satisfy 
condition  (11)),  the  auxiliary  output  cr(.r)  must  have  unitary 
relative  degree,  i.e.  its  first  derivative  must  explicitly  depend 
on  u.  Additionally,  the  signals  u+  and  u~  must  be  selected 
properly. 

Note:  In  the  case  of  linear  systems,  the  relative  degree  is 
given  by  the  difference  between  the  numbers  of  poles  and 
zeros  of  the  transfer  function. 

3.2.  Proposed  conditioning  technique 

3.2.1.  The  admissible  power 

The  specifications  of  the  electrolyzer  regarding  voltage 
limit,  current  limits  and  power  rate  limits  translate  into  lower 
and  upper  power  limit  signals  P,\(t)  and  P,\(t).  On  the  one 
hand,  the  upper  limit  signal  PaU)  is  determined  by  either  of 
the  rated  values  7n  or  Vn,  or  the  upper  power  rate  limit.  On  the 
other  hand,  the  lower  limit  signal  PaU)  is  established  by  /m 
or  the  lower  power  rate  limit.  It  is  defined  here  the  admissible 
power  PaU).  which  is  constructed  by  passing  the  electrolyzer 
current  through  a  limiter  upper-bounded  at  the  rated  current 
7n  and  lower-bounded  at  the  minimum  current  7m,  the  volt¬ 
age  through  a  limiter  upper-bounded  at  the  rated  value  Vn 
and  their  product  through  a  rate  limiter.  By  construction,  this 
admissible  power  ranges  over  PaU)  <  PaU)  <  PaU)  at  any 
time  t.  So,  the  specifications  for  the  electrolyzer  will  be  effec¬ 
tively  satisfied  whenever  the  actual  electrolyzer  power  equals 
this  admissible  power  (PuU)  =  Pa (!))■ 

On  the  basis  of  this  analysis,  the  following  output  signal 
is  defined: 


=  Ph  -  PaU).  (12) 

Note  that  whereas  (10)  defines  a  surface  in  conventional  slid¬ 
ing  mode  control,  in  the  current  application  o  =  0  deter¬ 
mines  a  region  in  the  state  space,  where  Ph  is  constrained 
within  its  admissible  limit  values  PaU)  and  PaU)-  in  fact, 
this  region  defined  by  TZ  =  {o  =  Ph  —  PaU)  =  0}  is  delim¬ 
ited  by  two  time-varying  surfaces  S  and  S  defined  by  A  = 
{o=  PH-  T\U)  =  0}  and  S  =  {a  =  PH  -  PaU)  =  0},  re¬ 
spectively. 

3.2.2.  The  SM  conditioning  objective 

Obviously,  the  specifications  of  the  electrolyzer  are  satis¬ 
fied  whenever  the  auxiliary  output  (12)  is  set  to  zero.  Thus, 
the  objective  of  the  conditioning  loop  is  to  maintain  the  sys¬ 
tem  operating  in  the  desired  region  TZ,  i.e.  to  constrain  Ph 
between  its  limit  values  PaU)  and  PaU)  all  the  time.  So, 
when  Ph  tries  to  cross  any  of  the  associated  limit  surfaces,  S 
or  S,  a  sliding  regime  is  established  until  Ph  is  able  to  evolve 
toward  the  interior  of  TZ  without  correction. 

In  contrast  with  conventional  SMC,  the  aim  of  the  con¬ 
trol  strategy  proposed  here  is  not  to  evolve  in  sliding  mode 


towards  the  equilibrium  point.  Contrarily,  the  sliding  regime 
is  intended  as  a  transitional  mode  of  operation.  It  is  aimed  at 
conditioning  the  power  reference  in  order  to  accomplish  the 
specifications  of  the  electrolyzer.  Hence,  once  the  captured 
power  evolves  smoothly  and  within  the  permissible  limits, 
the  sliding  mode  correction  becomes  inactive. 

3.2.3.  The  conditioning  loop 

The  conditioning  system  within  dashed  lines  in  Fig.  5  com¬ 
prises  a  block  devoted  to  compute  the  admissible  power  Pa 
as  detailed  above,  a  comparator  that  determines  the  auxiliary 
output  er,  a  switching  logic  that  decides  the  control  signal  u 
and  a  first-order  low-pass  filter  F  interposed  in  the  path  of 
the  power  reference  signal.  Note  that,  from  the  point  of  view 
of  implementation,  this  conditioning  system  can  be  realized 
as  a  set  of  instructions  at  the  software  level  of  the  MPPT 
controller.  So,  the  only  additional  hardware  is  the  on-line 
measurement  of  the  electrical  variables  of  the  electrolyzer. 

The  MPPT  reference  P*(T2)  is  affected  by  the  discontin¬ 
uous  factor  u  before  passing  through  the  low-pass  filter.  It 
must  be  remarked  that  this  filter  is  not  intended  to  smooth  by 
itself  the  reference.  In  fact,  this  obvious  approach  could  lead 
to  instability  as  explained  before.  Conversely,  the  proposed 
conditioning  lies  in  shaping  the  filter  output  Pr  as  function 
of  the  electrolyzer  specifications.  So,  the  cut-off  frequency 
of  the  filter  is  designed  to  attenuate  the  switching  frequency 
but  fast  enough  to  avoid  introducing  appreciable  lags.  The 
power  reference  Pr  is  thus  a  smooth  signal  shaped  by  the 
SM  control  loop  through  the  discontinuous  input  signal  u. 


3.2.4.  The  switching  logic 

Before  proposing  a  switching  law,  it  will  be  checked  the 
relative  degree  of  the  output  o  with  respect  to  the  discontin¬ 
uous  control  signal  u.  Note  that  o  and  u  are  the  output  and 
input,  respectively,  of  a  cascade  system  comprising  the  filter 
F,  the  PI  controller  and  the  ac-dc  converter.  The  filter  and 
PI  controller  are  linear  dynamical  systems  described  by  the 
transfer  functions 


F  :  F)s)  = 


1 

1  +  srF 


PI  :  C(s)  =  KP 


(13) 


where  rF  is  the  time  constant  of  the  filter,  K\>  the  proportional 
gain,  rc  the  integral  time  of  the  controller  and  s  is  the  Laplace 
operator.  Finally,  the  ac-dc  converter  may  be  described  by  the 
static  input-output  Eq.  (5).  Whereas  the  filter  F  has  unitary 
relative  degree  (it  has  one  pole  and  no  zero),  the  PI  controller 
(one  pole  and  one  zero)  and  the  converter  (its  outputs  V  j  [  and 
7h  explicitly  depend  on  its  input  S)  have  null  relative  degree. 
So,  the  cascade  system  F-PI—  ~  /  =  effectively  has  unitary 
relative  degree.  This  means  that  the  necessary  condition  for 
SM  existence  on  surfaces  S  and  S  is  satisfied.  Then,  sliding 
regimes  on  the  borders  of  the  desired  operating  region  TZ 
can  be  established  provided  an  adequate  switching  logic  is 
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decided  [16].  We  propose  here: 

{u  ~  =  1 ,  if  a  <  0 

u°  =  k,  if  a  =  0  (14) 

u+  =  0,  if  a  >  0 

Naturally,  the  sustainable  limits  for  the  conditioning  action 
are  optimum  power  (u  =  1)  and  no  power  (u  =  0).  The 
signals  u~  and  u+  are  designed  in  consequence.  Finally, 
u°  =  k  <  1  is  designed  as  a  compromise  between  efficiency 
and  available  control  effort.  In  fact,  with  regards  to  energy 
production,  the  parameter  k  should  be  selected  close  to  1 
so  that  Pr  ~  I1*  (52)  during  desired  operation  mode  (a  =  0). 
Nevertheless  a  factor  k  close  to  1  would  minimize  the  capa¬ 
bility  of  conditioning  the  electrolyzer  power  during  sudden 
wind  speed  decrease.  It  must  be  remarked  that  this  type  of 
compromise  commonly  appears  in  control  of  wind  energy 
systems  when  the  quality  of  the  power  supplied  to  the  grid 
or  load  is  an  important  issue.  In  our  application,  sacrificing 
some  percent  of  the  wind  energy  may  be  largely  rewarded 
by  the  improvement  of  the  electrolyzer  power  quality  as  well 
as  of  the  produced  hydrogen  one.  Actually,  since  the  turbine 
is  often  oversized  to  avoid  multiple  connections  and  discon¬ 
nections  of  the  electrolyzer  to  the  turbine,  some  reduction  of 
the  captured  energy  is  not  a  serious  disadvantage.  Moreover, 
the  available  wind  power  may  well  exceed  the  rated  power 
of  the  electrolyzer  during  long  periods  and  the  power  capture 
should  be  deliberately  further  reduced  for  safety  operation  of 
the  electrolyzer. 

3.2.5.  SM  conditioning  loop  operation 

The  conditioning  system  operates  as  follows.  Let  suppose 
that,  initially,  the  input  power  to  the  electrolyzer  satisfies  the 
specifications  and,  therefore,  u  —  u°  =  k.  Let  suppose  now 
that  the  wind  speed  (and  consequently  the  reference  power 
Pr)  increases  in  such  a  way  that  the  rated  values  of  cur¬ 
rent/voltage  at  electrolyzer  terminals,  or  its  power  rate  up¬ 
per  limit  is  reached  at  time  t.  At  that  moment,  the  actual 
electrolyzer  power  P\\(t)  reaches  its  maximum  permissible 
value  Ph(t).  Immediately  after  P\\(t)  tries  to  exceed  P,\(t),  u 
switches  to  u+  =  0,  thus  reducing  the  reference  power.  When 
this  occurs,  u  returns  to  its  default  value  m°  =  k,  so  P\\  will 
try  to  exceed  P,\(t)  again.  Consequently,  a  series  of  high  fre¬ 
quency  commutations  between  u  —  k  and  u  =  0  takes  place, 
establishing  a  sliding  mode  on  S  until  the  limitations  be¬ 
come  inactive  (Ph  =  P,\(t)  with  u  —  u°).  That  is,  when  Ph 
evolves  to  the  interior  of  the  desired  operating  region.  This 
occurs  when  the  voltage,  current  and  power  rate  fall  below 
their  rated  values  after  the  wind  speed  slows  down. 

Conversely,  when  Ph  tries  to  fall  below  Pa  (that  is,  when 
the  minimum  current  is  reached  or  the  captured  power  drops 
suddenly  because  of  a  decreasing  wind  gust),  the  factor  u 
switches  to  u~  =  1 .  In  fact,  the  controller  tries  to  increase  the 
supplied  power  in  order  to  restore  immediately  the  desired 
operating  mode.  When  this  is  achieved,  u  switches  again  to 
u°  =  k.  Asa  result,  a  sliding  regime  is  established  on  S,  with 


u  switching  between  u  =  k  and  u  =  1 .  When  the  limitation 
disappears  because  a  favorable  change  in  the  wind  conditions, 
the  normal  operation  (u  =  k)  is  recovered  and  the  electrolyzer 
variables  evolve  to  the  interior  of  their  permissible  ranges  of 
variation  without  sliding  mode  conditioning. 

Note  that,  independently  of  the  control  strategy,  the 
control  effort  that  is  available  to  attenuate  the  effects  of  wind 
drops  in  wind  energy  systems  is  limited.  In  fact,  the  energy 
stored  in  the  rotor  that  can  be  used  to  supply  the  power 
shortage  is  bounded.  In  the  current  application,  the  bounded 
control  action  means  that  the  sliding  mode  conditioning 
cannot  be  sustained  under  drastic  wind  reductions.  In  fact,  if 
the  wind  decreases  too  fast,  the  energy  stored  in  the  rotor  may 
not  be  sufficient  to  maintain  Ph  above  its  lower  power  rate 
limit  and  u  will  be  fixed  at  u~  =  1  trying  to  provide  as  much 
energy  as  possible  to  the  electrolyzer  during  the  transient. 
Similarly,  if  the  wind  decreases  to  a  too  small  value,  the 
available  wind  power  will  inevitably  be  not  sufficient  to 
supply  the  current  7m  to  the  electrolyzer.  This  limitation  is 
not  attributable  to  the  control  strategy  but  to  the  resource 
availability.  In  fact,  if  the  input  power  required  to  supply 
the  current  Im  to  the  electrolyzer  exceeds  the  available  in 
the  wind  (P0pt(uO),  maintaining  the  minimum  electrolyzer 
current  will  not  be  achievable,  independently  of  the  control 
strategy.  In  this  case,  the  backup  power  supply  should  be 
connected. 


4.  Simulation  results 

In  this  section,  simulation  results  are  presented.  They 
were  obtained  for  a  system  like  that  shown  in  Fig.  1  with 
a  2.25  kW  electrolyzer  and  a  5  kW  wind  turbine  plus  per¬ 
manent  magnet  induction  generator.  For  the  sake  of  clarity 
in  the  presentation  of  the  simulation  results,  the  electrolyte 
temperature  T  was  considered  constant.  Under  this  assump¬ 
tion,  current  and  voltage  bounds  lead  to  constant  power  limit 
signals. 

The  wind  speed  used  in  simulations  is  plotted  in  Fig.  6. 
This  wind  profile  has  been  intentionally  selected  to  evaluate 
the  performance  of  the  control  strategy  under  extreme  condi¬ 
tions.  The  results  of  three  simulation  runs  are  presented.  The 
first  two  sets  of  data  (cases  A  and  B )  were  obtained  without 
the  conditioning  loop,  whereas  sliding  mode  conditioning  is 
included  in  the  last  simulation  run  (case  C). 


Fig.  6.  Wind  speed  profile. 
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Fig.  7.  Simulation  results  with  MPPT:  (a)  electrolyzer  power  and  (b)  actual 
(thick  trace)  and  optimum  (noisy  trace,  overlapped)  turbine  speed. 

4.1.  Case  A:  maximum  power  tracking  control 

This  case  is  devoted  to  show  the  performance  of  the  con¬ 
ventional  MPPT  algorithm  for  comparison  with  the  proposed 
control  (case  C).  The  results  are  presented  in  Fig.  7.  Fig.  7a 
depicts  the  power  captured  by  the  turbine  and  supplied  to 
the  electrolyzer.  Fig.  7b  displays  the  optimum  speed  profile 
overlapped  by  the  actual  turbine  speed.  In  fact,  the  optimum 
speed  (for  maximum  power  capture)  is  successfully  tracked 
with  the  mechanical  dynamics  that  is  too  fast  to  be  distin¬ 
guished  in  the  time  scale  of  the  figure.  Thus,  the  electrolyzer 
power  also  overlaps  the  optimum  turbine  power  P0pt(w).  Note 
that,  although  the  optimum  operating  point  of  the  turbine  is 
successfully  tracked,  the  MPPT  performance  does  not  ac¬ 
complish  the  power  requirements  of  the  electrolyzer.  In  fact, 
the  supplied  power  reproduces  the  fast  wind  variations  and 
exceeds  during  some  time  intervals  its  rated  value  as  well  as 
its  power  rate  limit. 

The  option  of  smoothing  the  power  reference  P*(Q)  with 
a  low-pass  filter  is  evaluated  here.  Reference  filtering  is  often 
used  in  many  control  problems  to  avoid  nonlinearities  or  un¬ 
desirable  responses.  However,  this  approach  does  not  provide 
a  satisfactory  solution  in  this  application.  Fig.  8  shows  the 
simulation  results  obtained  when  the  power  reference  P*(P2) 
is  passed  through  a  first-order  filter  to  smooth  fluctuations 
and  then  limited  at  rated  power.  Results  are  presented  for 
two  cases:  Tc\  =  2  s  and  T(- i  =  3  s,  with  Tq  being  the  time 
constant  of  the  filter.  In  neither  of  the  two  cases,  the  filter 
smoothes  substantially  the  power  fluctuations.  Moreover,  in 
the  latter  case,  where  the  cut-off  frequency  is  lower,  the  filter 
introduces  a  lag  that  leads  to  a  turbine  stop.  This  undesirable 
behavior  has  been  predicted  in  Section  2.2.  Because  of  the 
filter  lag,  the  filtered  reference  often  exceeds  the  captured 
power  when  wind  drops,  hence  decelerating  the  rotor  below 
the  optimum  speed.  In  particular,  the  wind  gust  at  t  =  120  s 
drives  the  turbine  into  the  stall  region  (see  in  Fig.  8b  that 
the  rotor  speed  falls  well  below  the  optimum  speed,  i.e. 
X  <5IC  A()pt),  where  the  system  presents  an  unstable  behavior. 
Since  wind  varies  randomly,  there  is  no  guarantee  of  stability 


Fig.  8.  Simulation  results  with  filtered  reference  (dot-dashed:  tq  =  2  s, 
solid:  re  =  3  s,  thin:  without  filter):  (a)  electrolyzer  power  and  (b)  turbine 
speed. 

when  a  low-pass  filter  is  used  to  smooth  the  reference 
power. 

4.2.  Case  B:  suboptimum  power  point  tracking 

This  control  strategy  is  a  prelude  of  the  SM  condition¬ 
ing  technique  evaluated  in  the  following  paragraph.  It  basi¬ 
cally  consists  of  the  conventional  MPPT  algorithm,  where 
the  reference  power  locus  is  obtained  by  multiplying  P*{P2) 
by  a  coefficient  k  <  1,  more  precisely  k  =  0.7.  That  is, 
Pr  =  0.7  P*(Q).  Note  that  this  selection  of  k  does  not  im¬ 
ply  a  30%  reduction  of  the  captured  power  with  respect  to 
the  optimum  case.  In  reality,  the  system  searches  a  new  op¬ 
erating  point  tjk  at  higher  speeds,  where  the  power  reference 
curve  kP*(P2)  crosses  the  actual  power-speed  locus  of  the  tur¬ 
bine  (see  Fig.  9).  At  this  new  point,  the  turbine  operates  with 
a  power  coefficient,  lower  than  Cpopt  of  course,  that  depends 
on  the  shape  of  its  Cp-A.  curve. 

Fig.  10  shows  the  simulation  results.  Fig.  10a  illustrates 
the  captured  power  which  also  violates  the  specifications  of 
the  electrolyzer.  In  fact,  during  the  simulation  run,  the  sup¬ 
plied  power  exceeds  its  rated  value  (2.25  kW),  falls  below  its 


Fig.  9.  Operating  points  for  different  power  reference  loci. 
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Fig.  10.  Simulation  results  with  suboptimum  power  tracking:  (a)  elec¬ 
trolyzer  power  and  (b)  turbine  speed  (solid  line)  and  turbine  speed  for  case 
A  (dashed  line). 

minimum  permissible  value  (0.35  kW  in  this  example)  and 
presents  fast  variations.  Fig.  10b  reveals  that  the  turbine  is  ef¬ 
fectively  operating  at  a  speed  higher  than  in  case  A,  thereby 
reducing  the  energy  capture  as  expected.  Actually,  after  com¬ 
paring  the  captured  power  in  this  case  and  in  the  previous 
optimum  case,  it  follows  that  the  energy  loss  is  around  7%. 


the  sliding  regime  between  t  =  40  and  54  s  limits  the  decreas¬ 
ing  rate  of  change  to  its  lowest  admissible  value.  Also,  the  SM 
conditioning  between  t  =  78  and  1 14  s  successfully  confines 
In  (and  hence  P]\)  to  its  rated  value.  Finally,  from  t  =  365  to 
382  s,  the  sliding  mode  increases  the  reference  power,  thus 
avoiding  In  to  fall  below  Im.  Besides,  during  other  periods, 
the  electrolyzer  power  evolves  smoothly  and  within  its  per¬ 
missible  values  without  SM  conditioning.  See,  for  instance, 
the  response  from  t  —  382  to  430  s. 

Fig.  lib  illustrates  how  the  turbine  adjusts  its  speed  in 
order  to  supply  an  appropriate  amount  of  power  to  the  elec¬ 
trolyzer.  Note  that  the  speed  is  above  the  optimum  turbine 
speed  most  of  the  time,  particularly  during  regulation  at  rated 
power  or  maximum  power  rate,  whereas  it  approaches  £2opt 
during  regulation  at  minimum  power  or  minimum  power  rate. 

Fig.  lid  confirms  that  Pn  is  maintained  within  the  desired 
region  of  operation  (er  =  0)  all  the  time  except  for  a  short 
period  around  t  =  120  s.  During  this  period,  characterized  by 
a  fast  decreasing  gust,  the  conditioning  loop  applies  the  max¬ 
imum  control  effort  to  increase  as  much  as  possible  the  cap¬ 
tured  power.  Thus,  during  this  period  the  controller  behaves 
as  an  MPPT  control.  This  kind  of  behavior  is  not  completely 
avoidable  because  the  available  control  effort  is  limited.  The 
compromise  in  the  design  of  k  between  energy  capture  and 


4.3.  Case  C:  suboptimum  power  point  tracking  with 
sliding  mode  conditioning  to  satisfy  the  electrolyzer 
power  requirements 

The  simulation  results  obtained  with  the  proposed  con¬ 
trol  algorithm  are  shown  in  Fig.  11.  In  Fig.  11a,  the  power 
supplied  to  the  electrolyzer  Pn,  which  tightly  follows  the 
conditioned  reference  If,  is  plotted.  Fig.  1  lb  exhibits  the 
turbine  speed  response.  Fig.  11c  displays  the  discontinuous 
conditioning  signal.  Finally,  Fig.  lid  depicts  the  auxiliary 
output  response. 

The  simulation  results  corroborate  the  effectiveness  of  the 
proposed  conditioning.  In  fact,  the  electrolyzer  current  has 
been  bounded  below  its  rated  value  and  the  captured  power 
has  been  smoothed,  thereby  improving  the  quality  of  the 
power  supplied  to  the  electrolyzer. 

As  clearly  observed  in  Fig.  1  lc,  the  simulation  run  alter¬ 
nates  intervals  of  SM  conditioning  with  periods  of  no  cor¬ 
rection.  In  fact,  the  SM  conditioning  periods  are  identified 
by  the  fast  switching  behavior  of  the  discontinuous  signal  u. 
Moreover,  sliding  regimes  occur  on  the  upper-limit  surface 
a  —  0  when  u  switches  between  u°  and  0,  whereas  they  are 
established  on  the  lower-limit  surface  o  =  0  when  u  switches 
between  u°  and  1.  On  the  other  hand,  the  conditioning  loop 
is  inactive  during  the  nonzero  time  intervals,  where  u  is  fixed 
at  m°. 

Fig.  1  la  shows  how  the  conditioning  loop  shapes  the  elec¬ 
trolyzer  power  with  the  aim  of  accomplishing  the  specifica¬ 
tions.  For  instance,  from  t  =  0  to  34  s,  the  sliding  regime 
is  able  to  constrain  Pn  to  its  maximum  permissible  rate  of 
change  despite  the  rapidly  increasing  wind  power.  Similarly, 


Fig.  1 1 .  Simulation  results  with  suboptimum  power  tracking  and  SM  con¬ 
ditioning:  (a)  electrolyzer  power;  (b)  turbine  speed  (solid  line)  and  turbine 
speed  for  case  A  (dashed  line);  (c)  SM  control  signal;  (d)  SM  auxiliary 
output. 
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power  quality  improvement  is  evident.  Although  the  sliding 
regime  cannot  be  sustained  during  this  short  period,  the 
overall  response  is  much  better  than  in  the  unconditioned 
case. 

5.  Conclusions 

The  power  control  of  a  wind-hydrogenenergy  system  has 
been  addressed.  The  proposed  controller  basically  combines 
a  maximum  power  point  tracking  algorithm  (such  those  con¬ 
ventionally  used  to  maximize  efficiency  in  wind  energy  sys¬ 
tems)  and  a  novel  auxiliary  loop  that  suits  the  captured  power 
to  the  requirements  of  the  electrolyzer.  Actually,  this  loop 
shapes  the  power  reference,  when  necessary,  guaranteing  safe 
operation  of  the  electrolyzer  and  improving  the  conditions 
for  hydrogen  production.  The  reference  conditioning  circuit 
is  designed  using  concepts  of  sliding  mode  control  theory. 
The  proposed  solution  is  very  simple  to  implement  and  tune, 
and  shows  a  very  fast  response.  Furthermore,  it  can  be  readily 
incorporated  to  pre-existing  maximum  power  point  tracking 
controllers.  The  ability  of  the  controller  to  shape  properly 
the  power  profile  has  been  corroborated  by  a  simulation 
example. 
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